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Abstract. Q-LET is a FORTRAN 77 code that enables a quick estimate of the 
gravitational lensing effects on a point- or an extended source. The user provided 
input consists of the redshifts, angular positions relative to the source, mass or velocity 
dispersion estimate and halo type for the lens galaxies. The considered halo types are 
the Navarro-Frenk- White and the Singular Isothermal Sphere. The code uses the so- 
called multiple lens-plane method to find the magnification and intrinsic shape of the 
source. This method takes into account the multiple deflections that may arise when 
several mass accumulations are situated at different redshifts close to the line-of-sight. 

The Q-LET code is applied to the recently discovered supernova, SN2003es, which 
is likely to be of Type la as its host galaxy is classified as an elliptical. We find that 
SN2003es is likely to have been significantly magnified by gravitational lensing and 
that this should be considered in high-z studies if this SN is to be used to determine 
the cosmological parameters. 

Q-LET was motivated by the supernova searches, where lensing can be a problem, 
but it can also be applied to any simple lens system where a quick estimate is wanted, 
e.g. the single lens case. 
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1. Introduction 

The recent advance in modern cosmology due to new and improved telescopes and 
techniques seems to overthrow the long-lasting paradigm or hope of an Einstein-de 
Sitter (EdS) universe. From the Cosmic Microwave Background (CMB) ^EHHj, along 
with studies of Large Scale Structure (LSS), [4 and the Lyman a forest Ej (Ly^) 
etc. it is inferred that the universe consists of roughly 30 % matter (Cold Dark Matter, 
CDM and ordinary matter) and 70 % Dark Energy. Furthermore, measurements of the 
apparent magnitude of Type la supernovae (SNe) have been used by two collaborations 
[3 IHI to estimate the values of the total matter energy density ^m, which includes 
CDM, and the density in Dark Energy, e.g. ^a. Type la SNe are believed to be so- 
called standard candles, i.e. they show a very small scatter in intrinsic luminosity after 
empirical corrections jHHHII, and by measuring the apparent magnitudes, the cosmology- 
dependent luminosity distance can be inferred. The supernovae look fainter than they 
would in an EdS universe favouring both a non-zero f^A and f^M 7^ 1- Combined with 
CMB data from Boomerang and MAXIMA pTj IT^ . they end up at approximately 
the same values as the above-mentioned CMB, LSS and Lya measurements. However, 
there is a number of potential systematic errors that can contaminate the measurements. 
First, inter- and/or intra-galactic dust may obscure the light from the SNe and make 
them look fainter. This has been investigated in e.g. IT ^ ITH]. It is also still unclear 
whether there is any evolution of the SN brightness with redshift. Many potential sources 
of evolution have been proposed, such as differing progenitor composition jlHl HH UH] 
or host galaxy morphology (and thereby redshift) dependence ^Hl just to name a few. 
Furthermore, if photons can oscillate into axions, the SNe will also look fainter due 
to the reduced number of received photons. This has been addressed e.g. in reference 
PU] . A fourth possible contaminant is magnification or de- magnification by gravitational 
lensing. Having a large sample of SNe this can be corrected for statistically jSIj, but 
when studying single SNe, the lensing effects need a careful treatment. This was the 
case with the farthest known SN so far, SN1997ff, at a redshift of 1.7. The vicinity 
of the line-of-sight to this SN was unusually dense in galaxies and modelling of the 
individual lens galaxies was required. This was done by |22l 1211 who conclude a 
large possible magnification depending on the lens masses and concentrations. Benitez 
et al. conclude a magnification of 0.34±0.12 mag. 

In this paper we present the well-known multiple lens-plane method applied 
to the calculation of the magnification of any point- or extended source. Section 
|2] describes this method in a general way. Furthermore, we model the lenses as 
spherically symmetric Navarro-Frenk- White (NEW) or Singular Isothermal Sphere 
(SIS) halos and the lensing properties of these are described in sections 12.21 and 12.31 
For more information on under which circumstances each halo type is appropriate, 
see e.g. reference {21]. The method and models are manifested in a FORTRAN 77 
program, Q-LET, obtainable from the author upon request or it can be downloaded 
at http://www.physto.se/~cg/qlet/qlet .htm. The program is described in the text 
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in section IHl In reference j2^, an early version of Q-LET was used to estimate the 
magnification of the abovementioned SN1997ff. Section HI describes the apphcation of 
the code to the recently discovered supernova SN2003es. 

2. Multiple light deflection 

This section introduces the multiple lens-plane method and gives the important 
equations needed in the study of multiple light deflection. For more details, see reference 
|27j . We also describe the lensing properties of the NFW and SIS halo profiles. 

2.1. The multiple lens-plane method 

We start by projecting the mass of each massive object (e.g. galaxy) onto a plane at the 
respective object redshift. In the v.th. plane, this gives a surface mass density Sj(^j), 
where is the impact parameter of the light-ray in each plane (see figure Q). The 
justification for this projection is that the light is deflected only in the very vicinity 
of the lens in most cases and then travels practically unimpeded until it reaches the 
observer or another lens, situated far away from the preceding one. 

2.1.1. The lens equation From figure ^ where most quantities are defined, we see that 
the position of the light-ray in each plane can be obtained recursively from the observed 
position through 

i=l 

where Dij = D{zi,Zj) is the angular diameter distance between redshifts Zi and Zj, 
Di = D{0,Zi) and Dig = D{zi, Zgonrce)- With lens-planes we find the source position 
as T7 = ^7v+i- convenience we use, in each plane, dimensionless quantities defined 
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where is an arbitrary scale length in the i:th plane which, if chosen as Di, makes Xj 
the angular impact parameter in each plane. Furthermore, the convergence is 



Scrit 

where (in geometrised units) 
V ■ 



crit 



(3) 



(4) 



is called the critical density due to its relation to the ability of a lens to produce multiple 
images. If we also define cti = ctiDis/Dg, the scaled deflection potential ipi^Xi), needed 
for the magnification and defined through oci — Wipi becomes 



V'i(xi) = - / Kj(x') In |xj - x'|d^x'. 

71" Jr2 



(5) 



In these units where y = xat+i is the source position, the lens equation is given by 

N 

y = xi -^aj(xi). (6) 



2.1.2. The magnification We now want to find the magnification as a function of image 
position. Let A(xi) denote the Jacobian matrix of the lens equation 

dy 



Furthermore we need 



(7) 



/ dag dag \ 
' dxn dxi2 ' 



y dai2 dai2 I 
\ dxii dxi2 / 



where 



l^i + In 1x2 

li2 Ki - In 



dxkdxi ' 



(8) 



(9) 



(10) 



Vx,V'i = 2«;i(xj), 7ii = (ijjin - -0122)72 and 7^2 = V'a2 = -0121- This gives A as 

where I is the 2x2 unit matrix. The Afs can be found by recursion and by noting that 
A, = I- 

j-i 

Aj^I-J2WAi- (11) 

i=l 
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Finally, the magnification is given by 
1 



/i 



detA 



(AiiA 



22 



12^21; 



(12) 



Negative values of /i indicates images of reversed parity relative to the unlensed image. 



2.2. Properties of the NFW halo profile 

The use of the Navarro-Frenk- White model calls for a brief description of its lensing 
properties. 

Before projection onto the lens plane, the NFW profile is given by 



Ps 



1 + 



2 ' 



(13) 



where ps and are scale parameters that can be determined once e.g. the mass is 
specified, r is the radial co-ordinate. The projected surface mass density will be 
circularly symmetric and we no longer need vector notation on the impact parameter 
but instead choose ^ as the radial co-ordinate with the origin at the lens centre. If we 
choose ^0 = fs (cf. equation (0)) and define Hs = Psrs/^ait, then 



where 



k{x) 



-fix) 



1 
1 



arctan 



farctanh 



and k(1) = 2k,s/3. The deflection angle is 



x-l 
x+1 

1-x 
1+x 



,X > 1 
X <1 



a{x) = Aks 



X 



(14) 



(15) 



(16) 



where 



In - + 
2 



2 



arctan 
arctanh 
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x-l 
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1-x 
1+x 



X > 1 
X <1 
X = 1 



(17) 



and the magnification will be given by 

a{x) 



/i(x) 



a(x) , , 
1 + - 2k{x) 



X / \ X 

A flaw of this model is that the total mass diverges when integrated out to inflnity. 
However, this is not too serious since the deflection angle only is sensitive to the 
mass inside the impact radius, and the magniflcation is sensitive to this mass and the 
convergence at this point (for circularly symmetric lenses). Thus all mass outside the 
impact radius is unimportant. 

A single NFW halo gives either one or three images, where the primary image has 
p> I, the secondary < and the tertiary > 0. 
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2. 3. Properties of the SIS halo 

The Singular Isothermal Sphere model, also used here, is based on the assumption 
that the dark matter in the halo behaves as particles in an ideal gas trapped in their 
gravitational potential. The gas is assumed to be in thermal equilibrium and the 
resulting density profile before projection (in geometrised units) is 

Psis(r) = (19) 
where ct„ is the line-of-sight velocity dispersion of the particles. If we choose 

eo = 47rcr.^^ (20) 

the equations for convergence, deflection angle, intrinsic source position and 
magnification become very simple for this model: 

k{x) = ^, (21) 
a{x) = ^ (22) 



and 



and 



\x\ 



y{x) =x--^ (23) 
\x\ 



Kx) = yf^. (24) 
\x\ — 1 



As with the NFW model, the SIS will give an infinite mass when integrated to infinity 
but the argument above regarding deflection angle and magnification also holds here. 
When considering a single SIS halo it can give either one or two images depending on 
the impact parameter. Primary images have n > 1 and secondary have n < 0. 



3. The Q-LET code 

The Quick Lensing Estimation Tool, Q-LET is a FORTRAN 77 code written in order to 
quickly be able to estimate the lensing effects on a point source or an extended object. 
The structure of the code is as follows 

• Read the input datafile and sort the planes 

• Check if parameter values are allowed 

• Present a choice between a point- or an extended source with elliptical image shape, 
alternatively a square grid (see section 

• Compute magnifications and positions of light rays 

• Present a choice of how to output the results 
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A used supplied datafile contains the cosmological parameters, the source redshift and 
position and the lenses' redshifts, positions, velocity dispersions or masses and finally 
their halo type. The mass should be M200, i-e. the mass within the radius (= r2oo) within 
which the average energy density is 200 times the critical density at the lens' redshift. 

A point source must be put at the origin. For the case of an elliptical image, the 
origin will be put at the source centre and the ellipse will be centered at user specified 
co-ordinates which should be taken as (0,0) if an elliptical image only is studied. The 
possibility of putting the ellipse off-origin is useful to study e.g. a SN that is offset from 
its host galaxy centre and the galaxy shape is to be investigated. The SN is put at the 
origin and the host galaxy central co-ordinates are given as input. If a square grid (see 
section is used, it will be centered at (0,0). The angular positions of the lensing 
galaxies should of course be assigned with respect to this origin and should be given 
in arcseconds. The velocity dispersions and masses should be given in km s~^ and 
respectively. See section 13.21 for an example of an input file. 

If only M200 is given for a SIS lens, this mass will be used to calculate the 
corresponding velocity dispersion since is the parameter used to determine the lensing 
strength for this halo type. If, on the other hand, only is given for a NFW lens, this 
velocity dispersion will be used to calculate the corresponding M200 for a SIS lens and 
assuming M200 is the same for the NFW halo:}:. This is done since the lensing properties 
of the NFW halo depends upon M200 in our calculations. If, for some reason, both the 
mass and velocity dispersion are given, the mass (velocity dispersion) will be ignored 
for the SIS (NFW) lens. 

All distances are computed using the filled beam approximation and it is also with 
respect to an homogeneous universe that the magnification is given. 

3.1. P arameter restrictions 

Due to computational complications (and perhaps lack of connection to reality) there 
are some approximate restrictions on the input parameters. However, these restrictions 
do not severely limit the usefulness of the code as can be seen in table ^ where the 
allowed parameter ranges are presented. Furthermore, some of them are not very strict 
and depends on the values of the other parameters making these ranges approximate 
only. Especially extremely low or high Hubble parameters (~ 0.1 or ~ 1.5 might lead 
to problems. 

3.2. Input 

The information needed by Q-LET to be able to compute the magnification and 
defiection of a light ray is the following: the values of the cosmological parameters, 
i.e. /i, f2ni and ^a, the source redshift and position, the lenses' redshifts and positions, 
their velocity dispersion or mass and their halo types (NFW or SIS). This information 

X See section |S1 for a discussion of the validity of this. 
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Source redshift 


< < 20 


Lens redshift 


< < min[2;s, 5] 


Lens mass [M©] 


10^ < M200 < 10^^ 


Lens velocity disp. [km s~^] 


5< < 5000 


Hubble parameter 


< /i< 2 


Mass energy density 


< l^m < 10 


Energy density in cosmological const. 


-10 < f^A < 10 


Combined restrictions 


Jim + JIa = 1 or 


when using NFW halo 


< 1 and Qa = 



Table 1. Approximate parameter restrictions in the code. The Hubble parameter and 
energy densities are given in units of 100 km s~^ Mpc~^ and critical energy density, 
Pcrit respectively. 



is collected in an input datafile which could have the following appearance: 





0.7 


0.3 


0.7 










1.7 














XXX 


0, 


.923 


-4.1 


3.01 


138 





nf w 


0, 


.655 


2.6 


-1.11 





lel2 


sis 


0, 


.322 


-3.32 


5.11 





9.87el0 


nfw 


1, 


.453 


-2.32 


-0.4 


180.0 





sis 


0, 


.011 


-2.34 


2.32 





6.5e9 


nfw 


0, 


.234 


-4.1 


3.01 


190 





sis 


0, 


.254 


-5.23 


-2.43 


100.3 





nfw 



The entries are: 



-^source -^source [ ] Z/source [ ] 



XXX 



-^-lensi -^lensi [ ] Z/lensi 

["] (7vi [km s"^] M2001 [M©] halo typci (NFW/SIS) 



(7v„ [km s-i] M2oo„ [Mq] halo type„ (NFW/SIS) 

For easier handling of the data you should put '0' in the fourth and fifth columns 

and e.g. in the sixth column of row 2. Remember to put the source at (0,0) and of 

course at the highest redshift. The following rows (up to 199 more rows) should contain 
data on the lensing halos according to the example above. 

The program starts by asking for this datafile. The default name is 'qdata. qlet'. 
In the next step, a question is posed on whether you want to study a point source, an 
elliptical image or a 20" x 20" square grid§. If you choose an elliptical image, you will 

§ Can be used to trace an arbitrary image shape back to the source-plane to see the intrinsic shape. 



^le 



^lens„ ["] yiens„ ["] 




Figure 2. Definitions of major, a, minor, b, axes, rotation angle, (p and offset x' and 
y'. Units: arcseconds (axes and offset) and degrees (rotation angle). 



be asked for another input file with major and minor axes, rotation angle and possible 
offset from the source (see figureEl). The default name of this input file is 'elldata. ell' 
and consists of one line only. An example input file is: 
0.8 0.3 72.0 -0.47 0.12 

The inputs are a 6 in arcseconds, ip in degrees and the offset x' and y' in arcseconds. 
3.3. Output 

If a point source is chosen, it is placed at (0,0) and traced back to the source plane and 
the magnification is computed. This is printed on the screen, and it is also possible to 
obtain the output (x-, y-position and magnification) in a small file. If a magnification 
yU < 1 is obtained, it is not a primary image that has been found, i.e. muliple imaging 
has occurred. In the case of an elliptical image, the user is asked to provide the number 
of the plane wished to be output in a file. You can choose any plane between 2 and the 
number of the source plane if you would like to study the intermediate changes of the 
ray bundle. Each row in the output file refers to a specific ray and could look like: 
8395. 0.2131891 0.8527565 -2.296262 0.633064 2.873963 
The entries are: 

No. of the ray xi ["] yi ["] Xi ["] yi ["] magnification, 
where subscript 1 refers to the image plane and i refers to the chosen plane. The file will 
contain a quite large amount of rays in order to see the shape better (of order 1-5x10^, 
depending on the size of the ellipse). To see the result, just plot the file with the desired 
columns. 

When a square grid is chosen, the output file will look exactly like in the ellipse- 



Q-LET — Application to SN2003es 



10 



case, the difference is that the grid is square and pre-defined to be 20" x 20". This grid 
can be used if the image(s) has (have) a shape not described by an elhpse. When trying 
to estimate the intrinsic shape of an object like this, you plot the 4:th and 5:th columns 
making cuts using the 2:nd and 3:rd columns as to only include the rays making up the 
image(s). If instead wishing to see e.g. what an intrinsically elliptical source looks like 
in the image plane, you make the cut using the 4:th and 5:th columns, even though the 
grid in the source plane is not evenly distributed it will be possible as long as it covers 
the area where the cut is made. 



4. Example: application to SN2003es 

SN2003es was found at a redshift of 0.968 in the Hubble Deep Field North (HDF-N) in 
April 2003 by the GOODS/HST Transient Search, as part of the HST Treasury program 
|29] . As the host galaxy seems to be an elliptical, it is most likely a Type la SN |30] . 
Since these are used as standard candles to determine the redshift-distance relation in 
the quest for the cosmological parameters, it is very important to study the potential 
magnification by gravitational lensing. There are several possible lens candidates close 
to the line-of-sight of SN2003es so a close inspection is necessary. 



^.i. Host and lensing galaxy data 

The data on the host and lens galaxies has been taken from three different references; 
|3H IS^ IMj. Since the redshift information differ significantly in the different catalogs, 
we used the principle to take the spectroscopic redshifts when available and the 
photometric redshifts otherwise. If two catalogs had only photometric redshifts, the 
most recent result was used. However, as a worst case scenario regarding the amount 
of magnification, we also computed the magnification including four galaxies from one 
of the older catalogues which, in the most recent catalogue, were determined to be 
background galaxies aX z > 2;source- Hereafter we call these GH after Gwyn and Hartwick 
of reference [33 . 

To estimate the velocity dispersion that is needed in the lensing calculation, we 
took the observed luminosities and best fit galaxy morphology from and performed 
cross-filter K-corrections to rest frame Bj using galaxy templates provided by the SNOG 
package f^. The Tully-Fisher (T-F) relation for spirals (Sxx) and Faber- Jackson (F- 
J) relation for ellipticals (Ell) and irregulars (Irr) were used to derive an approximate 
velocity dispersion. The relations are given by 

^ _ ^q0.088(M.-M) _ (25) 

and 

^ = IQ^MM.-M) (p_j)^ (26) 

where cr* is the normalisation of the velocity dispersion, M is the absolute magnitude 
in Bj and M* is a typical magnitude taken to be M* = —19.53 + b\ogh 
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Cat. no 


Z X- 


-pos. ["] 


y-pos. ["] 




Gal. type 


526 


0.120 


-3.74 


-8.60 


0.136 


Scd 


539 


0.120 


-2.55 


-7.93 


0.112 


Scd 


549 


0.952" 


1.62 


-8.29 


0.924 


Ell 


622 


0.440 


4.76 


-2.59 


0.332 


Scd 


623 


0.960 


-6.80 


2.58 


0.334 


Scd 


631 


0.321^^ 


1.69 


-0.52 


0.460 


irr 




n '^1 1 " 


9 71 


9 RR 
z.oo 


(1 ^"^7 


Irr 




n finn 




u.oo 


(1 4(19 


£_/ll 


7nn 

( uu 


n (14(1 


1 ^9 




(1 (1Q8 


DDL 






TT 

GH g 


;alaxies 






582 


0.904 


-9.73 


-0.62 


0.401 


Scd 


607 


0.776 


1.64 


-3.07 


0.431 


Sbc 


610 


0.296 


0.87 


-2.43 


0.218 


Scd 


696 


0.904 


8.52 


4.51 


0.682 


Irr 






Host 


galaxy 






619 


0.968^^ 


-0.475 


-0.710 


(1.290) 


Ell 



Table 2. Lensing galaxy data. The galactic centre position is given relative to the 
co-ordinates of SN2003es. The galaxy type is the best fit from reference The 
galaxies with spectroscopically determined redshifts are superscripted a. 



Within a radius of 10" we found 9 lensing galaxies excluding GH and 13 lenses 
including it. The redshifts, velocity dispersions, positional data and galaxy types can 
be found in table |21 

4-2. Results 

We used the Q-LET code to compute the magnification of SN2003es and also to find the 
intrinsic shape of the host galaxy. This is done in order to see that the observed elliptical 
shape is not contrived from an intrinsically very complicated shape that is never observed 
in reality. The normalisation of the velocity dispersion is quite poorly known [HHI and 
we choose to present the magnification as a function of this normalisation. For the host 
galaxy shape calculation we used a value of a^, = 170 km/s, obtained as an average in 
reference |37j for a mixture of galaxy morphologies. However, the shape is not distorted 
very much even with cx^, as high as 300 km/s, although it is translated and more strongly 

II Catalogue number of reference 
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Figure 3. Magnification of SN2003es as a function of the velocity dispersion 
normalisation. 

magnified. The result on the magnification is found in figure El We have used models 
where all halos are assumed to be either of NFW or SIS type. The figure also shows the 
result both when including the four GH galaxies and when excluding them. As is seen, 
in the region of cr* between 150 and 200 km/s, i.e. around 170 km/s, the magnification 
varies between roughly 1.08 and 1.27 depending on the the model and the presence of 
the GH galaxies. This corresponds to a magnitude shift of 0.1 < Am < 0.25 mag, a 
significant change that could dominate the uncertainty (the intrinsic magnitude spread 
in Type la SNe is roughly between 0.1 and 0.2 mag [SH]). Note also that higher a* have 
been obtained e.g. in reference [39j, where a value for the SIS of roughly 275 km/s was 
obtained for the dark matter in ellipticals. 

In figure EJ we display the intrinsic (unlensed) position relative to the observed 
position at (0,0), depending on a^, and the model (NFW/SIS and GH incl./GH excl.). 
We see that the shift is only weakly depending on the halo type for this lens system 
under our assumptions. 

Finally, we show the intrinsic shape of the host galaxy for a a,, of 170 km/s jSZ] 
in figure El This plot also shows the small difference in deflection between NFW and 
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Figure 4. The intrinsic (unlensed) position of SN2003es relative to the observed 
position at (0,0). The position will move along the curves from left to right with 
increasing tr* as indicated. 



SIS halo models under our assumptions. The shapes were computed including the GH 
galaxies. The axis ratios of both the NFW and SIS case are roughly 0.9, making them 
perfectly normal El galaxies so no constraints on the magnification or model can be 
inferred from the host galaxy shape. 

5. Discussion — Q-LET & SN2003es 

One of the approximations that may be criticised is the fact that we extend the halos out 
to infinity. However, the lensing effect at large impact radii is very small. For SN2003es, 
the lenses with large impact radii contribute very little to the total magnification which 
is dominated by galaxies no. 631 and 662. 

Putting M200 of the NFW equal to that of the SIS for a given velocity dispersion 
is not obviously a good approximation. The most likely measured quantities are the 
luminosity or velocity dispersion of the lens and the luminosity may be used to estimate 
a velocity dispersion either via the TuUy-Fisher or Faber- Jackson relation. When a 
velocity dispersion is obtained, the best way would be to use it directly on the NFW 
profile. However, as opposed to the SIS, the NFW velocity dispersion depends on the 
radius. By following Lokas and Mamon in reference [2H] it is possible to obtain the 
line-of-sight velocity dispersion, crios, for well resolved sources, or, for small or distant 
sources, the "aperture velocity dispersion", (Tap, which is the average over an aperture 
centered on the object. We calculated both these quantities as a function of radius and 
found that Uap did not fall off dramatically. As an example, in a case where we put 



Q-LET — Application to SN2003es 



14 



^ 0.25 

u 
u 

^ 



-0.25 
-0.5 
-0.75 
-1 
-1.25 



Obs4i\»(l 



-1.5 




C. Gunnarsson. 2003 



NFW 



-1.25 -1 -0.75 -0.5 -0.25 0.25 0.5 0.75 

X [arcsec] 



Figure 5. The observed and intrinsic shapes of the host galaxy. The larger dots 
indicate the position of the supernova. The figure shows the small difference in 
deflection between NFW and SIS halo models under our assumptions. 

the SIS velocity dispersion (cr^^^) to 170 km/s to get the parameters for the NFW and 
then compute the corresponding NFW velocity dispersions, dap varies roughly by +5% 
and —15% around 170 km/s between radii 0.1r2oo and 2r2oo- The line-of-sight velocity 
dispersion varies more, +10% to —40% in the same range. However, at a radius of 
roughly 0.3r2oo, there is a quite small scatter around the SIS value of less than ~15%. 
This seems to be valid for a large range of redshifts and velocity dispersions. It also 
seems valid in low Hq- and low density universes. Thus in almost all cases it is possible 
to find a radius at which a\os ~ cxap ~ cTy^^) with this radius being roughly equal to 
0.3r2oo- 

Inverting the velocity dispersion-mass relation for the NFW is quite time- 
consuming, and this, along with the ambiguity due to the dependence of on radius 
and the fact that the difference seems not to be significant, motivates using the simpler 
conversion via the SIS. However, if this conversion is unsatisfactory, it is always possible 
for the user to supply M200 instead, in which case of course no conversion is needed. 

For the irregular galaxies we have used the Faber- Jackson relation which is valid 
for ellipticals. An approximation which is difficult to say to what extent it affects the 
results, and in which direction. Using the TuUy-Fisher relation on the irregulars would 
increase the magnification even further, but at most by ~ 15-20%. 

Whether the halo models we use are good descriptions of real halos is still 
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under debate. We are using spherically symmetric models, which might be an 
oversimplification but the ellipticity of the dark matter in the observed halos is in 
practice very seldom known (see e.g. reference [10] for the effects of ellipticity in lensing). 

So, unless some new model is very different from the SIS and NFW, the results for 
SN2003es is likely to remain. 

5.1. Conclusions 

We have found that SN2003es is likely to have been significantly magnified by 
gravitational lensing, a fact that must be taken into account when trying to use this 
(likely Type la) SN to infer cosmological parameters. Either a more thorough analysis of 
the lensing galaxies^ is needed to take the effect into account or the SN should be used 
only to set a lower limit on the distance or maybe even be excluded from the sample 
if not having large statistics. In that case some SNe are likely to be de-magnified and 
the effect can be assumed to average to ~ 0. The fact that the intrinsic shape of the 
host galaxy looks perfectly normal makes it practically impossible to use this shape to 
put constraints on the magnification of SN2003es. If an extremely unusual shape would 
have been obtained, it could perhaps have affected the reliability of the model. 

6. Summary 

We have presented the multiple lens-plane method of gravitational lensing along with 
its implementation in Q-LET, a FORTRAN 77 code that can be used to quickly estimate 
the lensing effects on either a point source or to obtain the intrinsic shape of an elliptical 
image. A square grid in the image plane is also possible for user-defined image shapes. 

The code was applied to SN2003es, a SN (most hkely Type la) discovered in the 
HDF-N. Using both NFW and SIS type halos and estimating the velocity dispersions 
from the observed luminosities via the Faber- Jackson and Tully-Fisher relations we 
found that SN2003es is likely to have been significantly magnified by gravitational 
lensing and that this problem should be addressed in detail before using SN2003es 
in estimates of the cosmological parameters. 
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